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We investigate the possibility to explain theoretically the galaxy rotation curves by a gravitational 
potential in total absence of dark matter. To this aim an analytic fourth-order theory of gravity, 
nonminimally coupled with a massive scalar field is considered. Specifically, the interaction term 
is given by an analytic function f{R,<j)) where R is the Ricci scalar and (p is the scalar field. The 
. . . , gravitational potential is generated by a point-like source and compared with the so called Sanders's 

CO ' potential that can be exactly reproduced in this case. This result means that the problem of dark 

' matter in spiral galaxies could be fully addressed by revising general relativity at galactic scales and 

, requiring further gravitational degrees of freedom instead of new material components that have not 

been found out up to now. 
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I> '• 

I ^! I. INTRODUCTION 

O ; 

^ '. J.H. Oort first pointed out the missing matter problem in the 30's of last century [l], [2[. The issue came out by 
^ ' observing the Doppler shift of stars moving near the plane of our Galaxy and calculating the star velocities. The 
result was that there had to be a large amount of matter inside the galaxy to prevent the stars from escaping. Such a 
"matter" should give rise to a central gravitational force much larger than Sun's gravitational pull to keep a planet in 
T-H its orbit. However it turned out that there was not enough luminous mass in the Galaxy to account for this dynamics. 
^ The discrepancy was very large and the Galaxy had to be at least two or three times more massive than the sum of all 
its luminous components in order to match the result. Later on, the tangential velocity of stars in orbits around the 
Galactic center was calculated as a function of distance from the center. Surprisingly it was found that far away from 
the Galactic Center, stars move with the same velocity independent of their distance out from the Galactic Center. 
• These results strongly posed the problem that either luminous matter was not able to reliably trace the radial profile 
ii^ ' of the Galaxy or the Newtonian potential was not able to describe dynamics far from the Galactic center. 

, Soon after, other dark matter issues came out from dynamical descriptions of self-gravitating astrophysical systems 
I ' like stellar clusters, galaxies, groups and clusters of galaxies. In all these cases, there is more matter dynamically 
• • ' inferred than that can be accounted for by luminous matter components. The mass discrepancy comes out assuming 
. !^ the validity of Newton law at any astrophysical scales. Problems emerged also at larger scales. F. Zwicky discovered 
' anomalous motions of galaxies in the Coma cluster finding that the visible mass was too little to produce enough 
■ gravitational force to hold the cluster together [Sj]. 
. At the beginning, the only possibility considered was to assume the Newton law holding at all scales and postulating 

some non-luminous component to make up the missing mass. Many names have been coined to define these invisible 
components. For example, the MAssive Compact Halo Objects (MACHOs) are objects like black holes and neutron 
stars (in general sub-luminous objects) that populate the outer reaches of galaxies like the Milky Way. There are the 
Weakly Interacting Massive Particles (WIMPs) which do not interact with standard matter (constituted by baryons 
as protons and neutrons): they are supposed to be particles out of the Standard Model of Particles but, up to now, 
there is no final indication for their existence In general, dark matter is assumed to come in two flavors, hot 
(HDM) and cold (CDM) dark matter. The CDM should be in dead stars, planets, brown dwarfs etc., while HDM 
should be constituted by fast moving relativistic particles. It should be neutrinos, tachyons etc. However, there is 
still no definitive proof that WIMPs exist, or that MACHOs will ever make up more than five percent of the total 
amount of missing matter. 

On the other hand, the need of unknown components as dark energy (coming from cosmology) and dark matter 
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could be considered nothing else but as a signal of the breakdown of Einstein General Relativity (GR) at astrophysical 
(galactic and extragalactic) and cosmological scales. 

In this context, Extended Theories of Gravity (ETGs) could be, in principle, an interesting alternative to explain 
cosmic acceleration and large scale structure without any dark components. In their simplest version, the Ricci 
curvature scalar R, linear in the Hilbert-Einstein action, could be replaced by a generic function ,f{R) whose true 
form could be "reconstructed" by the data. In fact, there is no a priori reason to consider the gravitational Lagrangian 
linear in the Ricci scalar while observations and experiments could contribute to define and constrain the "true" theory 
of gravity (see [Bl-flHV 

Coming to the weak-field limit, any alternative relativistic theory of gravity is expected to reproduce GR results 
which, in any case, are firmly tested only at Solar System scales in the Newtonian limit Even this limit is 

matter of debate since several relativistic theories do not reproduce it. For example, Yukawa-like corrections to the 
Newtonian potential easily comes out [isj with interesting physical consequences. For example, it is claimed by some 
authors that the flat rotation curves of galaxies can be explained by such terms [TJ]. Other authors have shown that 
a conformal theory of gravity is nothing else but a fourth order theory containing such terms in the Newtonian limit. 
In general, any relativistic theory of gravitation yields corrections to the weak- field gravitational potentials {e.g., [Tsl l) 
which, at the post- Newtonian level and in the Parametrized Post-Newtonian formalism, could constitute a test of 
these theories 

This point deserves a deep discussion. Beside the fundamental physics motivations coming from Quantum Gravity 
and unification theories (see [1, iHl] ) , ETGs pose the problem that there are further gravitational degrees of freedom 
(related to higher order terms, non-minimal couplings and scalar fields in the field equations) and gravitational 
interaction is not invariant at any scale. This means that, besides the Schwarzschild radius, other characteristic 
gravitational scales could come out from dynamics. Such scales, in the weak field approximation, should be responsible 
of characteristic lengths of astrophysical structures that should result confined in this way [l6j . 

In this paper, without claiming for completeness, we will try to address the problem of describing galaxy rotation 
curves without dark matter but asking for corrections to the Newtonian potential that could fit data and reproduce 
dynamics. These corrections are not phenomenological but come out from the weak field limit of general relativistic 
theories of gravity that predict the existence of corrections (e.g. Yukawa-like corrections) to the Newtonian potential. 
The only exception is GR where the action is chosen to be i?, that is linear in the Ricci curvature scalar and does not 
contain corrections to the Newtonian potential in the weak field limit. Relaxing such a hy pothesis, it is possible to 
show that any analytic ETG presents Yukawa corrections in the weak-field limit (see also |15| for a detailed caculation). 
From an astrophysical point of view, these corrections means that further scales have to be taken into account and 
that their effects could be irrelevant at local scales as Solar System. With this scheme in mind, we will give a summary 
of ETGs in Sec. |ll] discussing also their conformal properties. In fact any ETG can be conformally transformed to 
the Einstein one plus scalar fields representing the further gravitational degrees of freedom. This feature is extremely 
important to select characteristic length scales (related to the effective masses of scalar fields) that could account for 
dynamics. In this sense, considering /(i?) gravity means to take into account an Einstein theory plus a scalar field; 
considering f{R, □i?)-gravity means to assumes Einstein -I- two scalar fields and so on. The emergence of Yukawa-like 
corrections to the Newtonian potential is discussed in Sec. IIIII where the weak-field limit of /(i?)-gravity, the simplest 
ETG, is worked out. Here, /(i?) is a generic analytic function of the Ricci curvature scalar R. Furthermore, we discuss 
the case of /(i?, 0)-gravity, corresponding to /(i?, □-R)-gravity, i.e. Einstein plus two scalar fields, showing that a 
further free parameter in needed to better model dynamics. Sec. IIVI is devoted to the rotation curves of galaxies. It 
is shown that the phenomenological Sanders potential, suitable to fit realistically observations, can be reproduced by 
the weak filed limit of /(i?, (/))-gravity. Sec|V]is devoted to discussion and conclusions. 



II. EXTENDED GRAVITY AND CONFORMAL TRASFORMATIONS 

Higher-order and scalar-tensor gravities are examples of ETGs. For a comprehensive discussion, see (5l-[To|. Essen- 
tially these theories can be characterized by two main feature: the geometry can non-minimally couple to some scalar 
field; derivatives of the metric components of order higher than second may appear. In the first case, we say that 
we have scalar-tensor gravity, and in the second case we have higher-order theories. Combinations of non-minimally 
coupled and higher order terms can also emerge in effective Lagrangians, producing mixed higher order/scalar-tensor 
gravity. The physical foundation of such models can be found at fundamental level by considering effective actions 
coming from quantum fields in curved space-times, string/M theory and so on ^T^. A general class of higher-order- 
scalar-tensor theories in four dimensions is given by the effective action 

5 = / d^x^ [fiR, OR, a'^R, a'^R, 0) + u;((^)(/.;„<^^" + XCm] , (1) 
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where / is an unspecified function of curvature invariants and scalar field (j) and X = SirG^. The convention for 
Ricci's tensor is i?^^ — R'^fj.au, while for the Riemann tensor is i?"/3^i/ = T^^^ + .... The affinities are the usual 

Christoffel symbols of the metric: F^^ = ^g^'^{gaa,i3 + gi3cr,a — 9ai3,cr)- The adopted signature is (H )). The 

term Cm is the minimally coupled ordinary matter contribution, considered as a perfect fluid; Lo{(j)) is a function of 
scalar field which specifies the theory. Actually its values can be lo{4)) = ±1,0 fixing the nature and the dynamics 
of the scalar field which can be a canonical scalar field, a phantom field or a field without dynamics (see jl7Hl9| for 
details). In the metric approach, the field equations are obtained by varying IT} with respect to f^^^. By introducing 
the Einstein tensor G^u we get 

i=l j=l ^ ' \^ ^ ' 

where we have introduced the quantity 



J=0 



the energy-momentum tensor of matter 



1 6{y/~^C n 

~g Sgf^" 



T,. = -^^^^^^T^ (4) 



and □ = ;cr''^ is the d'Alembert operator. The differential Eqs.(l2]) are of order at most (2fc + 4). The (possible) 
contribution of a self- interaction potential V{(t)) is contained in the definition of /. By varying with respect to the 
scalar field (/>, we obtain the generalized Klein-Gordon equation 

2uj{^)U^ + u^icp) - /0 = , (5) 

where ftp = and uj^ {(j)) — ^^^^^ . Several interesting cases can be worked out starting from the action ([l} . Below, 
we give some significant examples that will result useful for the astrophysical applications of this paper. 

A. The case of /(i?)-gravity 

The simplest extension of GR is achieved assuming. 



and the action ([T} becomes 



Then the field equations ^ become 



R^fiR), L.(^) = (6) 



S = J d^xV^ [f{R) + XCm] ■ (7) 



/fl Gpi. = XTf^v + - — Y^di^y + fR:j_iu - g^yOfR = XT/j,^ + fnT^l^'' (8) 



^ Here we use the convention c = 1. 
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df . . . . . 

where fu = — — . The gravitational contribution due to higher-order terms can be reinterpreted as a stress-energy 
uR 

tensor contribution T/^i, . This means that additional and higher-order terms in the gravitational action act, in 

principle, as a stress-energy tensor, related to the form of /. In the case of GR, r/i^'' identically vanishes while the 
standard, minimal coupling is recovered for the matter contribution. 

The peculiar behavior of f{R) ~ R is due to the particular form of the Lagrangian itself which, even though it 
is a second order Lagrangian, can be non-covariantly rewritten as the sum of a first order Lagrangian plus a pure 
divergence term. The Hilbert-Einstein Lagrangian can be in fact recast as follows: 



-g = 



(9) 



where: 



n"/3 — 



-99 



dC 

dRaf 



(10) 



r is the Levi-Civita connection of g and u'^^ is a quantity constructed out with the variation of T [20|. Since uj^^ 
is not a tensor, the above expression is not covariant; however standard procedures can be used to recast covariance 
[lOj . This clearly shows that the field equations has to be of second order and the Hilbert-Einstein Lagrangian is thus 
degenerate. 



B. The case of scalar-tensor gravity 

From the action ([l}, it is possible to obtain another interesting case by choosing 



/ = F{^)R + V{(l>) , u;{^) = 1/2 , 



(11) 



then 



S ^ 



F{(I>)R + ¥{(()) + 



(12) 



where V{(^) and F(0) are generic functions describing respectively the potential and the coupling of a scalar field 
The Brans-Dicke theory of gravity is a particular case of the action (fT^ for 1^(0) = [l^. The variation with respect 
to g^j^y gives now the second-order field equations (particular form of field equations ([5])) 



F(0) G^, = XT^y + ^ff^. + F{4>),^, - 9t..^F{cj,) - i ( cj>,^cj>,y 



, ^ XT,, + F{<t>)Tl£> (13) 



where T^i^"* is the energy-momentum tensor relative to the scalar field 0. The variation with respect to (j) provides the 
Klein - Gordon equation, i. e. the field equation for the scalar field 



where F^(0) 



dF{cj>) 



dV{cj,) 



U(j)-F^{(j))R-V^{(l))^Q, 



. This last equation is equivalent to the Bianchi contracted identity [21| 



(14) 



C. Conformal transformations 



These models, and, in general, any theory of the class ([T|), can be conformally reduced to the Einstein theory 
plus scalar fields. Conformal transformations are mathematical tools very useful in ETGs in order to disentangle the 
further gravitational degrees of freedom coming from general actions [22| - [2^ . The idea is to perform a conformal 
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rescaling of the space-time metric g^^^ g^i^. Often a scalar field is present in the theory and the metric rescaling is 

accompanied by a (nonlinear) redefinition of this field — > 0. New dynamical variables (f^ are thus obtained. 

The scalar field redefinition serves the purpose of casting the kinetic energy density of this field in a canonical form. 

The new set of variables l^^^jy, ^| is called the Einstein conformal frame, while {gfj,^, 4>} constitute the Jordan frame. 

When a scalar degree of freedom (/) is present in the theory, as in scalar tensor or f{R) gravity, it generates the 
transformation to the Einstein frame in the sense that the rescaling is completely determined by a function of 0. In 
principle, infinitely many conformal frames could be introduced, giving rise to as many representations of the theory. 

Let the pair {M,gp,u} be a space-time, with Ai a smooth manifold of dimension n > 2 and g^^ a (pseudo)- 
Riemannian metric on M . The point-dependent rescaling of the metric tensor 

ffAii. — > 9tiv = ^^gtiu , (15) 

where the conformal factor £7 is a nowhere vanishing, regular function, is called a Weyl or conformal transformation. 
Due to this metric rescaling, the lengths of space-like and time-like intervals and the norms of space-like and time-like 
vectors are changed, while null vectors and null intervals of the metric g^i, remain null in the rescaled metric 5^^. 
The light cones are left unchanged by the transformation ([T5|) and the space-times {TW,^^^} and {A^,^^^} exhibit 
the same causal structure; the converse is also true [1^. A vector that is time- like, space-like, or null with respect to 
the metric 17^1, has the same character with respect to g^^, and vice-versa. 

Conformal invariance corresponds to the absence of a characteristic length (or mass) scale in the physics. In general, 
the effective potential V{(l)) coming from conformal transformations contains dimensional parameters (such as a mass 
m, that is a further "characteristic gravitational length"). This means that the further degrees of freedom coming 
from ETGs give rise to features that could play a fundamental role in the dynamics of astrophysical structures. In 
what follows, we will see that these further gravitational lengths could solve, in principle, the dark matter problem. 

D. Conformal transformations and higher-order gravity 

Performing the conformal transformation for /(i?)-gravity with — fji we have 

J d^x^g[f{R) + XLra\ = y d^x^g {r + W{4>) ~ + Xt^ (16) 

where = \fi In //j while the potential W and the nonniinimally coupled lagrangian of ordinary matter are 
given by 

(17) 

The function V is defined by the analogy between the /(i?)-gravity and the scalar-tensor gravity (the so-called 
O'Hanlon lagrangian) 

V(c^) ^ f[R)-RfR{R) (18) 
where (j) — fj^. The field equations in standard form are given in the Einstein frame as follow 

G^u = X H ^ — g^,y + - I (p-^j,(p-u — 9py I (19) 

n~^ + w.C^) = -x^-^ (20) 
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However, the problem is completely solved if ((> = Jr can be analytically inverted. In summary, a fourth-order theory 
is conformally equivalent to the standard second-order Einstein theory plus a scalar field (see also [26l. ^t\). 
If the theory is higher than fourth order, we have Lagrangian densities of the form |28l - l30l |. 

C = C{R, OR, n'^R). (21) 
Every □ operator introduces two further terms of derivation into the field equations. For example a theory like 

C = RDR, (22) 
is a sixth-order theory, and the above approach can be pursued considering a conformal factor of the form 

" dR^ dDR- ^^'^> 

In general, increasing two orders of derivation in the field equations [i.e. every term Di?), corresponds to add a scalar 
field in the conformally transformed frame j29j . A sixth-order theory can be reduced to an Einstein theory with two 
minimally coupled scalar fields; a 2n-order theory can be, in principle, reduced to an Einstein theory + (n— l)-scalar 
fields. On the other hand, these considerations can be directly generalized to higher-order-scalar-tensor theories in 
any number of dimensions as shown in (lo| . With these considerations in mind, we can easily say that a higher order 
theory like f{R, Di?) is dynamically equivalent to f{R, 0). This feature, as we will show, gives the minimal ingredients 
to reproduce the rotation curves of galaxies since two Yukawa like corrections come out. For a detailed derivation see 

m. 



III. YUKAWA-LIKE CORRECTIONS TO THE GRAVITATIONAL POTENTIAL 



In order to deal with standard self-gravitating systems, any theory of gravity has to be developed to its Newtonian 
or p ost-Newtonian limit depending on the order of approximation of the theory in terms of power of velocity 
|2lL l3]|. The paradigm of the Newtonian limit starts from the development of the metric tensor (and of all additional 
quantities in the theory) with respect to the dimensionless velocity^ v of the moving massive bodies embedded in 
the gravitational potential. The perturbative development takes only first term of 0,0- and i, j-component of metric 
tensor 17^1, (for details, see [32l. Issj). The metric assumes the following form 



ds^ = (1 + 2$) dt^ - (1 - 2*) S.jdx'dx^ (24) 

where the gravitational potentials $, vj/ < 1 are proportional to w^. The adopted set of coordinates^, the so-called 
isotropic coordinates, is x'' = (i,x) = {t,x^ ,x^,x^). The Ricci scalar is approximated as i? = i?'-'^' -1- i?^^) -I- . . . 
where R^^^ is proportional to $, and 4', while R'^^^ is proportional to 4'^ and 

Here we show as a general gravitational potential, with a Yukawa correction, can be obtained in the Newtonian 
limit of any analytic /(i?)-gravity model. From a phenomenological point of view, this correction allows to consider 
as viable this kind of models even at small distances, provided that the Yukawa correction turns out to be not relevant 
in this approximation as in the so called " chameleon mechanism" (33 | . 



A. Yukawa-like corrections in /(7?)-gravity 



Starting from the action (P) for the case f{R, DR, D^i?, . . . , D'^R, (/))-|-w(0)(/).c0'" reduced to f{R), the field equations 
are ([8|). In principle, the following analysis can be developed for any ETGs. Let us now start with the f{R) case. As 
discussed in IsH |36| . we can deal with the Newtonian limit of f{R) gravity adopting the spherical symmetry. By 
introducing the radial coordinate r — |x| the metric (l24l) can be recast as follows 



^ The velocity v is expressed in unit of light speed. 

^ The Greek index runs from to 3; the Latin index runs from 1 to 3. 
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ds^ = [1 + gil^ {t, r)] dt^ - [1 - (t, r)] dr^ - dQ , (25) 

where dil — dO^ + sin'^ 9 dcjP' is the angular distance, (i, r, ^, </>) are standard coordinates. Since we want to obtain 
the most general result, we do not provide any specific form for the f{R). We assume, however, analytic Taylor 
expandable f{R) functions with respect to the value i? = (Minkowskian background): 



h + fiR 



n = 



6 



(26) 



In order to obtain the weak field approximation, one has to insert expansions (|25p and (|26|) into field Eqs. ([8]) and 
expand the system up to the orders 0(0) and 0(1). This approach provides general results and specific (analytic) 
theories are selected by the coefiicients fi in Ea.p6|. It is worth noticing that, at the order 0(0), the field equations 
give the condition fo — Q and then the solutions at further orders do not depend on this parameter as we will show 
below. If we now consider the 0(1) - order approximation, the field equations in vacuum (T^j/ = 0), results to be 



Tglt 



(1) 



2/2ri?(i) = , 



(1) 



0, 



(27) 



/iri?(i) + 3/2 2i?W + ri?(i 



25. 



(1) 



= 0. 



It is evident that the trace equation (the fourth in the system (I27p ). provides a differential equation with respect to 
the Ricci scalar which allows to solve exactly the system ([27)) at 0(1) - order. Finally, one gets the general solution: 



(1) _ . Y 5,{t)e 

gtt - <^o - 7 — I — z — 5~ 

jir dm'^r 



_^ (52(t)e" 



Y 6i{t)[l + mr]e- 



fir 

^(1) ^ Si{t)e- 



+ 



Smrr 
2mr 



62{t)[l-mr\ 



(28) 



where = —^j^, and Y are arbitrary constants, while 5i{t) and 52{t) are arbitrary time functions. When we 
consider the limit f{R) — J> R then m — ?> 00 and /i — > 1, in the case of a point-like source of mass M, we recover 
the standard Schwarzschild solution if we set 60 — and Y — 2GM. Let us notice that the integration constant Sq 
is dimensionless, while the two functions 6i{t) and S2{t) have respectively the dimensions of length~^ and length~'^. 
These functions are completely arbitrary since the differential equation system (j27p contains only spatial derivatives 
and can be fixed to constant values. Besides, the condition 5q = is avaible since it represents an unessential additive 
quantity for the potential. 

The solutions ([28]) are valid if > i.e. fi and /2 are assumed to have different signs in Eq. ([26]) . If the 
algebraic signs are equal we find an oscillating solution where the correction term to the newtonian component 
(oc 1/r) is proportional to the (cosmr -I- smmr)/r (32j . In this paper we consider only the correction Yukawa-like. 

It is possible now, to write the general solution of the problem considering the previous expressions (|28p . In order 
to match at infinity the Minkowskian prescription for the metric, one can discard the Yukawa growing mode in 
and then we obtain : 
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1 - 



2GM ^ Si{t)e- 



R 



8i{t)e- 



2GM Siit)[l + mr]e- 
fir im?r 



dr^ 



r^dQ . 



(29) 



At this point, one can provide the solution in terms of gravitational potentials. The first of ([28)) gives the first order 
solution in term of the metric expansion (see the definition (j25p V This term coincides with the gravitational potential 

at the Newton order. In particular, since gtt — I + 2^grav = 1 + gtt \ the gravitational potential of /(i?)-gravity, 
analytic in the Ricci scalar _R, is 



fGM 5i(i)e-™'-\ , , 

V fir Gm'^r ) 

This general result means that the standard Newton potential is achieved only in the particular case /(-R) = R while 
it is not so for analytic /(i?) models up to exceptions of measure zero. Specifically all models with /1//2 > are 
excluded by hand. Eq. (j30l) deserves some comments. The parameters /i, m and the function (5i(t) represent the 
deviations with respect the standard Newton potential. To test these theories of gravity inside the Solar System, we 
need to compare such quantities with respect to the current experiments, or, in other words. Solar System constraints 
should be evaded fixing such parameters (33 |. On the other hand, these parameters could acquire non-trivial values 
(e.g. /i ^ 1, (5i(i) ^ 0, m < 00) at scales different from the Solar System ones. Since the parameter m can be related 
to an effective length i^^, Eq. ([50)1 can be recast as 

*-- = -iTT^j r ' (3^) 

M 

where the first term is the Newtonian-like part of the potential for a point-like mass _^ ^ ^ and the second term is 

a modification of gravity including a new scale length, L associated to the coefficients of the Taylor expansion. If 
(5 = the Newtonian potential and the standard gravitational coupling are recovered. Comparing Eqs. ([50)1 and 

pip , we assumed /i = 1 + 5 and (5i(t) = — ( r ) where 6 can be chosen quasi-constant. Under this 

L-^ \1 + J 

assumption, the scale length L could naturally arise and reproduce several phenomena that range from Solar System 
to large scale structure [16|. Understanding on which scales the modifications to GR are working or what is the 
weight of corrections to Newton potential is a crucial point that could confirm or rule out these extended approaches 
to gravitational interaction. 



B. Yukawa-like corrections in /(i?, (^)-gravity 

A further step is to analyze the Newtonian limit starting from the action ([ij and considering a generic function of 
Ricci scalar and scalar field. Then the action becomes 



A = 



(32) 



The field equations are obtained from ^ by setting /(i?, Di?, D^i?, . . . , □'^i?, 0) f{R, (f>)- As discussed in Sec. II, 
this case can be considered from a pure ly g eometric point of view assuming f{R,DR) theories where the terms OR 
gives a further scalar field contribution fl5|. We get 



JrR^u 7, 9^J.l' + ^{<P) 9;fi 4>;u - jB-f.i' + gt^u^-i JR = X 1 



2 w(0)n0 -I- tj^(0) -/^ = 



(33) 
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A further equation is the trace of field equation with respect to the metrci tensor g^,y 

fn R - 2/ - uj{<j)) + 3nfR = XT (34) 

where T — T" ^ is the trace of energy-momentum tensor. 

Let us consider a point-like source with mass M . The energy-momentum tensor is 

Tp^^ = pu^u^, T = p (35) 

where p is the mass density and satisfies the condition g^^uo^ — 1 and Ui — 0. Here, we are not interested to the 
internal structure. It is possible to analyze the problem in the more general case by using the isotropic coordinates 
{t,x^,x'^,x^), then the metric is expressed as in Eq. (j24p . In this framework, also the scalar field (j) is approximated 
as the Ricci scalar. In particular we get (p — 0^°) -f- (f)^^'> + (j)^"^^ + . . . and the function /(i?, (j)) with its partial 
derivatives (/^j, Jrr, f^, f^^ anf /^r) and Lu{(j)) can be substituted by their corresponding Taylor series. In the case 
of /(i?, (f>), we have 

f{R, cf>) ^ /(O, <^(o)) + fR{0, 0(°)) + UiO, <^(o)) 0(1) + . . . (36) 
and analogous relations for the derivatives are obtained. From the lowest order of field Eqs. (|33p we have 



/(O,0W) = O, /^(O,0W)-O (37) 

and also in this modified fourth order gravity a missing cosmological component in the action (1) implies that the 
space-time is asymptotically Minkowskian (the same outcome of previous section); moreover the ground value of scalar 
field (j) must be a stationary point of potential. In the Newtonian limit, we have 



A 



A 



/fli^(O,0^°)) 
/fl,(0>W) 

/fl(O,0(°)) 



Ac 



Ai?(i 



/fi0(O,0^°)) 
/i?(O,0(o)) 



i?(i) 



,fe4O,0(°)) ( 
/i?(O,0(°)) ^ 



1) 



j?(i) _ 
2 

i?(i) 



Xp 



/i?,(0, </'(")) 



/«(O,0(«)) /i.(O,0("))^ 



2t^(0(")) 

/fl(O,0(°))J?(i) ^ 
3/flfl(O,0(o)) 



i?(i 



/i?.0(O,'/'^°^) 
2tj((/.(o)) 

_ fR.Ao,<l>^'^) . , (1 

3/flfl(O,0(o)) /flfl(O,0(o)) 



= 



(38) 



where A is the Laplacian in the flat space. The last equation in ([55]) is the trace coming from Eq. ([M)) . These 
equations are not simply the merging of field equations of /(i?)-gravity and a further massive scalar field, but are due 
to the fact that the model f{R, (j>) generates a coupled system of equations with respect to Ricci scalar R and scalar 
field (j). By supposing that /^^ ^ and obviosuly /rr ^ 0, we can introduce the two characteristic length scales 



rriR 



3/flfl(O,0(o))' 



/00(O,'^'"^) 



(39) 



where the two masses are assumed to be real and this gives further restrictions to the set of viable models. The 
gravitational potentials $ and are given by*^ 



4 The potential "I- can be found also as = ^ f d^x'^^^^ + 4^ " •''"""^^"'tm' ' ^(^'(x). 
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ci>(x) 



X 



4^/i?(O,0(o)) 



1 

8^ 



' |x-x'| 3m«2 + /;,X0,<^("') ^ ^ 



(40) 



while, for the Ricci scalar and the scalar field, we have the coupled systera of equations 



A - 



A - 



(1) _ _MM!^i?(i) 
^ " 2^^(0(0)) 

(1) ^ mfl^A-p 3mHVfl^(O,0W) 
/h(O,0(°))^ /i?(0,<^(«)) 



(41) 



The definition of mji' is the generalization of w? in the case of pure /(i?)-gravity. By using the Fourier transformation 
the system (|4T|) has the following solutions 



p(k) 



i?(i)(x) = 



2c.(</)(0)) /fl(0, 0(0)) 7 (27r)3/2 (k^ + fci2)(k2 _^ ^^2) 



mj^^A" /■ d3k p(k) (k^ + gik-: 



(42) 



/fl(0, 0(0)) 7 (2^)3/2 + k^){\.^ + fcs') 



where 



2 2^1 3/^40, 0(")) m^^ ^ 



2L.(0("))/fl (0,0(0)) 



2^ 2 3/fl0(O,0(o))^m^2 



2^(0(o))/fl(O,0(o)) 



(43) 



In order to understand the relevant physical consequences of solutions (PO)) . it is sufhcient to analyze the point-like 
source framework. Then, if we consider p(k) — M j ^ where M is the mass, and fci ^'^ > Eqs. become 



0(1) (x) = 



/fl^ (0,0(0)) 



iRfcl |X| 



lRfc2 |x| 



2^.(0(0)) /fl,(0, 0(0)) 



^2 „ 1,2 



(44) 



i?(i)(x) = - 



"T-i? "Tg (kf - 77 )e 



2^ p-mjjfci |x| _ ^^2 _ 7y2^ g-mj7fc2 |x| 



/«(O,0("))|x| 



i-2 1,2 
'^1 '^2 



where is the Schwarzschild radius. Furthermore, we introduced the dimensionless quantities 



2^1,2 ~ 



2fcl,2 

m_R2 



= 1 - e + ± v/(i - e + ^')' - V 



(45) 



and 77 = — —, t, 



3/^^^^(0,0(0))^ 



mfl'^ 2a;(0(o))/^(O,0(o)) 



These two parameters have to ensures two different conditions for the roots 



fci_2 that have to be both real and positive i.e. fci.2 > 0. Such conditions can be reformulated as ^ < (77 — 1) and 
this fact restrict the class of viable Lagrangians. In fact we have 



tur'^ 



(46) 
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where ^ and rj are given in terms of fci.2 and mu which are the parameters defining the form of Yukawa-hke terms 
in the potentials. Specifically the conditions 



(47) 



/^40, = A/-C.(</>("))/fl(0>(0)) 



(fci^ - mfl2)(fc2^ - mij2) 



which, together with conditions ([37l) . give the form of possible Lagrangians. It is worth noticing that fii{0,4>^^^) can 
be assumed equal to 1 in standard units, while fixes the form of scalar field kinetic term that is equal to 1/2 

in the canonical case. For aj((/)'^*'^) < a ghost scalar field is possible. 
The potentials pO|) become 



$(x) 



*(x) 



GM 



/k(0>W)|x| 
GM 



where giCiV) 



/fl(O,0(o))|x| 
fc2(2772-2fc2-2^ + 3)-3r?2 



1 + V) e-"-^^ I''! + [1/3 - v)] e"™"^^ I"' 



(48) 



3^1 (^1 ^i) 

of the scalar field, the Ricci scalar and the potentials $, 



In Figs. ([T]), (l2|), (|3]) and (|4]) we show respectively the spatial behaviors 



4 6 
ijibIxI 



FIG. 1: The spatial behavior of scalar field <j!>'^-' generated by a point-like source (|44p for rj = 0.1 and ^ = —2. 

The solutions (|44| . (|48l) are the generalization of solutions obtained in /(i?)-gravity and scalar tensor gravity. In 
fact, we can easily obtain the outcomes that in the case of minimally coupled scalar field i.e. f^^ = 0—5-^ = 0, 
= 1; 9{S.tV) = 1/3 we find the point-like solutions of /(i?)-gravity [3^.[33j 



/«(0)|x|l ^3 



./«(0)|x|l 3 



(49) 



while in the case of the Brans-Dicke theory i.e. futf, = 1, /_r = 4>, (jj{4') = — wo/0, /rr = 0, /^^ = 0, from the 
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«Jb|x| 



FIG. 2: The spatial behavior of scalar field R^^^ (dashed line) generated by a point-like source (144 |l compared with respect to 
the same quantity (dotted line) in the /(7?)-gravity. In both cases we set rj = 0.1 and ^ — —2. 



-0.4 



O -0.8 



mslxl 



FIG. 3; Comparison among the potentials $ generated by a point-like source in three frameworks: the first one (|48|) induced 
by the action H32[) (dashed line), the second one (dotted line) induced by /(i?)-gravity (|49|) and the last one (solid line) is the 
Newtonian limit of GR. In the two alternative theories, we set rj = 0.1 and ^ = —2. 



field Eqs. we find the classical solutions of Brans-Dicke gravity |37| 

GM 2(2 + Wo) 



*i3D(x) = - 



(/.(o)|x| 2 Wo + 3 

GM 2(1 + Wo) 
0(0) |x| 2 Wo + 3 

In fact we recover the solutions (jSOp when mn 



(50) 



The Brans-Dicke behavior is also present in the solutions 
3 

oo, = 0, ^ = . Also solutions (j49|) have been found in vacuum, but here also the boundary conditions for 

2wo _ ^ 

the Yukawa term in the origin (where the mass is placed) have been inserted. In this case, the arbitrary time-function 

2m2 



6i{t) in Eq. (pO| is fixed to the value 
in the Eq. ([2^ becomes 



/i 



and then we can define $ 



grav 



while the expression for 



grav 



yrr 



GM 
1^ 



1 - 



(51) 



It is possible to show that the potential ^grav is equal to the potential ^ of Eqs. if we assume standard 

coordinates. The passage from the isotropic coordinates {t,x^,x'^,x^) to the standard ones {t,r,9,(j)) is given by the 
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FIG. 4: Comparison among the potentials ^ generated by a point-like source in three frameworks. The scheme is the same of 
Fig. El 



transformation 



l-2vl//(«)(|x|) 



(or rg/|x|), the metrics 



r- where I then, at first order with respect to the quantity rg/r 



//?(0)|x| 



'g 
hr 



-ni J?|x| 



dt^ - 



dt^ 



1 



g-™n|x| 



Sijdx^dx-' 



hr 



1 - 



1 -|- mr 



dr^ - r^dfl 



(52) 



coincide and, obviously, for f{R,(j)) — >■ /(i?), we have mn = m and /fl(0) = /i. 

It is interesting to note that, in the case of minimally coupled scalar field {fB.4> — 0), the Newtonian level of the 
field Eqs. p3p for the metric tensor is unaffected by the presence of the scalar field (/). Moreover is not linked to 
the energy-momentum tensor via the Ricci scalar and must satisfy only the boundary conditions at infinity, while the 
amplitude of scalar field is generic and depending only on time. In fact, the solution of the first of Eqs. (|4T|) is 



0(1) (t,x) 



K{t) 



-m^|x| 



2w((/)(o)) 



(53) 



where K{t) is a generic function depending on time. The evolution of K{t) is fixed by the post-Newtonian level of 
field equations. By considering //j^ ^ 0, we find a further contribution in the energy-momentum tensor. Another 
interesting case is the generalization of Brans-Dicke theory. In fact, we can consider a scalar-tensor theory, but the 
geometric sector is given only by the Ricci scalar. Without losing generality, we can set the interaction term in the 
action (j32p as (f) R and not as /(0) R because by introducing a new scalar field we obtain formally the same equations 
[11]. By setting /^^ = 1, Jrr = 0, Jr = cj) the field Eqs. §Eil become^ 



^ Also in this case we can find the solution of ^ as Nl/fx) = J- f d^x'-? — ^-^-X — ^A^^ 
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^ $ 



A - TO^^ 



(54) 



0(1) = 



2w(0(o)) 

r(i) - 3A0(i) 
0(0) 0(0) 



and their solutions are 



^ 2w(0(o))0(o)-3|x|^ 



i?(i)(x) = _1^5(x) 



6^(0(0)) .-V ,:riro;:iC. -^i-i 

[2w(0(o))0(o) -3]2 |x| 



(55) 



V 2 .,(0(0)) 0(0) -3 *' ' 

0(o)|x|\^ 2 a; (0(0)) 0(0) -3 

0|x| 



^St(x) 



GM r e V 2^(0(o))0{o)_3 



0(o)|x|\ 2a;(0(o))0(o) -3 



which, in the case of massless scalar field, become the typical solution of Brans-Dicke theory. 

In the cases that we have shown, the Newtonian contribution (|x|^i) to the potential is ever present. We can 
find a difference in the definition of gravitational constant G, since in these theories we have a multiplying factor 
/fl(0, 0(°))~i, while the additional terms are depending on the form of the Lagrangian. It is important to stress that, 
in all cases that we have considered, the limit and results of GR are fully recovered. 

If we have a generic matter source distribution p(x), it is sufficient to use the superposition principle by starting 
from point-like solutions. Then we substitute to the solutions the integral expression: $ — J This approach 
is correct only in the Newtonian limit since such a limit correspond also to the linearized version of the theory. 



IV. ROTATION CURVES OF GALAXIES 



At astrophysical level, the probe for the validity of alternative theories of gravity is the correct reproduction of 
rotation curves of spiral galaxies [1^1 ■ As discussed above, the foundation of the dark matter issue lies on this 
observational evidence. In order to face such a problem, one has to discuss the motion of a body embedded in a 
gravitational field. Let us take into account the geodesic equation 



where ds = gapdx°-dx^ is the relativistic distance. In the Newtonian limit, from Eq. (j56p . we obtain the equation 
of motion equation 



(57) 
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In our case, the gravitational potentials are given by (|48)) . The study of motion is very simple if we consider a 
particular symmetry of mass distribution p, otherwise the analytical solutions are not available. Our aim is to 
evaluate the corrections to the classical motion in the easiest situation: the circular motion. In this case we do not 
consider the radial and vertical motions. The condition of stationary motion on the circular orbit is 



«c(|x| 



9£(x) 
' a|x| 



(58) 



where Vc is the velocity. Generally the correction terms do not satisfy the Gauss theorem f39l| and this aspect implies 
that a sphere cannot be simply reduced to a point. In fact the gravitational potential generated by a sphere (also 
with constant density) is depending also on the Fourier transform of the sphere [s^. Only in the limit case, where 
the radius of the sphere is small with respect to the distance (point-like source), we obtain the simple expression (|48p . 

A further remark on Eqs. (|48|) is needed. The structure of solutions is mathematically similar to the one of fourth- 
order gravity generated by /(i?, -Ra/ji?"'^). However there is a fundamental difference with the present case: the two 
Yukawa corrections have different algebraic sign. In particular, the Yukawa correction induced by a generic function 
of Ricci scalar implies a stronger attractive gravitational force, while the second one, induced by squared Ricci tensor, 
implies a repulsive force [32|. |40j. In the present paper, the Yukawa corrections are induced by a generic function of 
Ricci scalar and a non-minimally coupled scalar field. Both corrections have a positive coefficient. In fact in Fig. ([5]), 
we show the coefficient g{S,,r]) with respect to for a given values of rj. The function g{^,r]) assumes the maximum 
value ( = 1/3) when ^ = (we have a pure /(i?)-gravity and the scalar field does not contribute to the gravitational 
potential in this case), otherwise we have two Yukawa corrections with positive coefficients. The scalar field gives rise 
to a more attractive force than in /(i?)-gravity. The interesting range of values of i] is between and 1. In the case 
7y > 1 — > 7710 > mji, the correction induced by scalar field is suppressed with respect to the other one. 




FIG. 5: Plot of coefficient g{^, rj) with respect to quantity ^ for < 77 < 0.99 with step 0.33. 



From this analysis, the choice of /(i?, 0)-gravity is better than /(i?, i?Q^i?"'^)-gravity, but we have a problem in 
the limit for |x| 00: the interaction is, of course, scale-depending (the scalar fields are massive) and in the vacuum 
the corrections turn off. For this reason, at large distances, we recover only the classical Newtonian contribution. 
Therefore the presence of scalar fields makes the profile smooth. This behavior is very clear in the study of rotation 
curves ((55)) . Let us assume a phenomenological point-like gravitational potential as supposed by Sanders [T^ 



$sp(x) = -^(1 -Ha e-™l-l) (59) 
|x| 

where a and 777 are free parameters that, following Sanders [jlj . can be assumed to be a ~ —0.92 and rp — 1/m ~ 
40 Kpc to fit the galactic rotation curves. This potential has been introduced to explain the rotation curves of spiral 
galaxies [4ll . |4^ . however the theoretical framework generating it is purely phenomenological. Recently by using 
the same potential it has been possible to fit elliptical galaxies [i^. In both cases by setting a negative value to a 
an almost constant profile of rotation curve is recovered. Such a rotation curve is obviously possible but there are 
two problems: the first one consists that no f{R, 0)-gravity, by imposing all boundary conditions at origin and at 
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infinity, gives that negative value of a. The second one is linked to the value of gravitational constant G. In fact 
in presence of Yukawa- like correction with negative coefficient, we find a lower rotation curve and only by resetting 
G (or the point-like mass) we can fit the experimental data. In Fig. Q we compare the profiles derived in the 
Newtonian limit of GR, f{R)- and /(i?, 0)-gravity and the potential (15^ . It is extremely interesting to note that 
the presence of scalar the field cf), in the case ~ m^, guarantees a rotation curve higher than the other ones but 
also in this case we find the asymptotic flatness are derived from observations. Only if we consider a massive scalar 
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FIG. 6: Comparison (in the vacuum case) of the spatial behaviors of rotation curves in the gravitational field generated by a 
point-like source. The dotted and dashed lines are the Sanders model for —0.95 < q < —0.92, the solid line is the GR curve, 
the dotted line is the /(i?)-gravity and the dashed line is the fourth order gravity non-minimally coupled with a scalar field. 
In the last case, we set ^ = —2, r\ = 0.1. 



tensor theory non-minimally coupled, we get a potential with negative coefficient in Eg. (1551) . In fact by setting the 

gravitational constant as Gn — "^/^x ^ where Gno is the gravitational constant as measured at infinity 

^ 2w(</>(")) -3 0(0) 

and by imposing a^^ = 3 — 2^(0'^°)) the potential $5t in the ([55]) becomes 



(60) 



and then the Sanders potential ([5^ is fully recovered. 



V. DISCUSSION AND CONCLUSIONS 



The dark matter issue, together with dark energy, can be considered the major problem of modern astrophysics 
and cosmology. Beside the huge amount of observations confirming its effects, practically at all astrophysical scales, 
no final answer exists, at fundamental level, definitively confirming one (or more than one) candidate supposed to 
explain the phenomenology. Furthermore, GR has been firmly tested only up to Solar System scale and then its 
features have been inferred at larger scales. In this situation, dark matter and dark energy could be nothing else but 
the manifestation that GR does not work at IR scales. 

A similar disturbing situation is found at UV scales where no Quantum Gravity theory is up to now definitely 
available. Alternative gravities (in particular ETGs) could represent a way out to this puzzle being effective theories 
of gravity representing a reliable picture of quantum fields in high curvature regimes [lo| and an approach to overcome 
the dark side problem at larger scales. 

In this paper, we have discussed the weak field limit (in particular the Newtonian limit) of some classes of ETGs in 
view to explain the almost fiat rotation curves of spiral galaxies. In particular, we have shown that ETGs, in general, 
present Yukawa- like corrections in the gravitational potential. In particular, we have analyzed the case of /(i?) and 
/(i?, 0). The latter are known to be analogue to /(i?, Di?). 

After a discussion of the mathematical features of the emerging corrections, we have confronted the results with 
the phenomenological Sanders potential, assumed as a possible dynamical explanation of flat rotation curves. The 
suitable value of phenomenological parameters can be exactly reproduced in the framework of /(-R, 0)-gravity since 
the concurring Yukawa corrections allow to recover attractive and repulsive components of potential. In this case, no 
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dark matter is required to fit dynamics like in the case discussed in (4^ where only a Yukawa-like correction was not 
sufficient to reproduce realistic rotation curves. 
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